Spin- Wave Relaxation in Diluted Magnetic Semiconductors within the Self-Consistent 

Green's Function Approach 
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We employ a self-consistent Green's function approach to investigate the spin-wave relaxation F(p) 
in diluted magnetic semiconductors. We find the trend of the spin-wave relaxation strongly depends 
on the ratio of the itinerant and impurity spin densities. For density ratios in the Ruderman-Kittel- 
Kasuya-Yosida phase, F(p) decreases even though thermal fluctuations increase. On the other hand, 
in the strong coupling phase, an interesting peak structure appears. We discuss the implications of 
our numerical results for experiments. 

PACS numbers: 75.40.Gb, 75.50.Dd, 76.50.+g 



Carrier-mediated ferromagnetism, found in diluted 
magnetic semiconductors (DMS) such as III-V semicon- 
ductors doped with transition metals, opens up the pos- 
sibility of electric manipulations of magnetic and optical 
properties because of the exchange coupling between the 
localized moments of the transition metal and the itin- 
erant spins in the semiconducting bands 0, 0] ■ In recent 
years, experimental progress has stimulated intense the- 
oretical investigations of DMS 0, S H H @ , which 
mainly focus on estimates of the critical temperature. 
Although some experimental data is available [3, E3 , lit- 
tle theoretical attention has been paid to spin dynamics, 
which is crucially important in understanding the spin 
coherence. 

To analyze how spin waves relax in DMS, both the itin- 
erant and impurity spins need to be treated on an equal 
footing. In this Letter, we employ the self-consistent 
Green's function approach[3 13 El [3 13 03 , which 
retains both spatial and thermal fluctuations, to inves- 
tigate the spin- wave relaxation at finite temperatures. 
The Green's functions allow us to numerically compute 
the spin spectral function and therefore the spin-wave 
relaxation. One may naively expect a smooth increase 
in the spin-wave relaxation as thermal fluctuations grow. 
However, our numerical results go against common in- 
tuition. For instance, in the Ruderman-Kittel-Kasuya- 
Yosida (RKKY) phase, the spin relaxation, shown in 
Fig. ^ actually decreases with increasing temperatures. 

To study the temperature dependence of the spin- wave 
relaxation, we begin with the simple Zener model, 




H = H + J J d 3 r S(r) ■ er(r), 



(1) 



where S(r) and <x(r) are the impurity and itinerant spin 
densities, coupled by the antiferromagnetic (J > 0) ex- 
change coupling. For simplicity, the kinetic energy of 
the itinerant carriers Hq is assumed to be parabolic, 
e p = p 2 /2m*. We take the typical values[l?J for the 
exchange coupling J = 150 meV nm 3 and the impurity 



FIG. 1: Spin spectral function A(p, Q,) / (2(S Z )) at tempera- 
tures T = 2,28K with density ratio nh/ni = 0.3. Against 
common intuition, the half-width shrinks as the temperature 
increases. 



spin density ni = 1 nm 3 but vary the ratio of itinerant 
to impurity spin densities n^/nj in the following. 

We now introduce the self-consistent Green's function 
approach. Note that, instead of treating spin waves as 
bosons, the Green's function approach respects the spin 
kinematics which is crucially important when thermal 
fluctuations are not small. The thermal Green's func- 
tion, describing the (impurity) spin-wave propagation, 
is defined as iD(r u r 2 ;t) = ((T[S + (ri,t)S-(r 2 ,0)])), 
where T is the time ordering operator and the double 
bracket implies both quantum and thermal averages. To 
complete the self-consistency, it is necessary to intro- 
duce another Green's function which describes the cor- 
relations between the impurity and the itinerant spins, 
iF{r u r'i,ri;t) = ((Tty}(ri, iJ^W, t)S~(r 2 , 0)] )). To 
simplify the calculations, we assume that the disorder is 
not strong, i.e. not close to the percolation threshold, and 
the virtual crystal approximation is valid. Therefore, af- 
ter a coarse-graining procedure, the translational invari- 
ance is approximately restored and the Green's functions 
only depend on the relative distance. This approximation 
is justified by recent numerical calculations |l8j which 
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FIG. 2: The magnetization curves of impurity and itinerant 
spins at different density ratios rih/ni = 0.01, 0.1, 0.3. 



found that magnetization curves with different disorder 
configurations are almost identical. 

The dynamical equations for the Green's functions can 
be derived by the standard approach and higher-order 
Green's functions may be decoupled within the random 
phase approximation. Making use of the translational in- 
variancc and the decoupling scheme, the dynamical equa- 
tion for D(r;t) in the Fourier space is 

nD( P ,n) = 2(s z ) - j(<r z )D(p,n) 

/H 3 k 
—F(k,k+p;n), (2) 

where F(k, k + p;£l) is the Fourier transformed Green's 
function F(r,r';t). Similarly, the other dynamical equa- 
tion is F(k, k+p;tt) = G(k, k + p; Q,)D(p, n) with 

G(k,k +P ;n)= J - /t(^-/i( e * + p) (3) 
2 S I + e k - e k+p + A + i-f h 

where f a (e k ) = [ e ^+ aA / 2 -^) _ jj-i ig thc Fcrmi digtri _ 
bution function for itinerant carriers with the Zeeman 
gap A = J(S Z ). On substituting the expression for 
F(k,k + p; fi) into Eq. (J2J, the spin-wave propagator is 
solved, 



DM = 



2{S Z 



(4) 



where the prefactor 2(S Z ) comes from the exact treat- 
ment of spin-wave kinematics. The self-energy E(p, f2) 
arises from interactions between the itinerant and the 
impurity spins and contains two terms, 



s(p,n) = -j{<t z 



d 3 k 



G(k,k+p;Cl). (5) 



If only the first term is retained, it coincides with the 
Weiss mean-field approximation. However, since the first 
term is always real, dropping the second term would re- 
move all information about the Landau damping of thc 
spin waves, described by the imaginary part of the self- 
energy Ej(p, f2), inside the Stoner continuum. 



Once the spin-wave dispersion is obtained from io v — 
E^(p, Up) = 0, the polarization of the impurity spins can 
be computed by Callen's formula [l9(, 

(S z )_ q _, (2S+l)(n sw y s +i 

\ n Bw/-r M , /„ \\2S4-1 /- \ 2.9+1 ' W 



m 



(1 + (n sw )) 2S + 1 - {n 8W ) 



where the average number of spin waves is {n sw ) = 
{1/ni) J d 3 fc/(27r) 3 [e^" - The difference between 

the independent spin- wave theory and the self-consistent 
Green's function method lies in the third term, which 
correctly accounts for the spin kinematics. The mag- 
netization curves for different density ratios rih/ni are 
shown in Fig. |5J While this is not the main focus of our 
results, we emphasize the sensitivity in the shape of the 
magnetization curve to the density ratio. From previous 
studies in the diffusive regime, there are three phases for 
the Zener model - mean-field, RKKY and strong cou- 
pling. For n^/nj = 0.3, it belongs to the RKKY phase 
and the magnetization curve is rather conventional. For 
nh/ni — 0.1, it lies in the strong coupling phase and the 
magnetization curve decreases linearly with a very small 
gradient and, only in the narrow regime near the criti- 
cal temperature, the magnetization dives to zero and the 
phase transition occurs. Similar magnetization curves, 
quite different to what we expect for a Heisenberg-likc 
model, have been observed in some experiments |20| . Fi- 
nally, we choose rih/nj = 0.01 (also in the strong coupling 
phase) to demonstrate that the magnetization curve can 
turn concave, though this is not necessarily related to 
percolation near the metal-insulator transition. 

The unique feature of the Green's function approach 
is that we can also study the temperature evolution of 
the spin spectral function, A{p,uj) = — (l/7r)ImD(p, O), 
as shown in Fig. ^ From Eq. Q , the imaginary part of 
the spin- wave propagator takes the Lorentzian form, 



A(p,u J ) = 2(S z ) 



nj (n-LU p ) 2 +T(pf 



(7) 



where Z w 1 is the spectral weight for the gapless spin 
wave. After Taylor expansion of E/(p, Q) in the vicinity 
of the spin- wave dispersion f2 = uj p , the relaxation rate 
of the spin waves T(p) w E/(a> p ). The numerical results 
for the spin- wave relaxation rate versus temperature with 
nh/ni — 0.3,0.1 are shown in Fig. EI One may expect 
that increasing thermal fluctuations lead to a increasing 
relaxation rate, but our numerical results disagree. In 
the RKKY phase with nh/ni = 0.3, the spin relaxation 
rate T(p) decreases as the temperature increases. This 
reduction can be understood as the weakening of the ef- 
fective coupling strength between the itinerant and im- 
purity spins in proportion to the spin polarization. How- 
ever, since the decoupling scheme in the self-consistent 
approach only includes decaying channels in the diago- 
nal parts (as in all random phase approximations), the 
computation of T(p) breaks down in the vicinity of the 
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FIG. 3: Spin-wave relaxation rate T(p) at different momenta 
and density ratios nh/ni = 0.3,0.1. The solid (broken) lines 
correspond to the left (right) axis. 

critical temperature where off-diagonal channels domi- 
nate. Therefore, one should not take the vanishing T(p) 
at the Curie temperature seriously. However, in the low 
temperature regime, the diagonal channels do dominate 
and our numerical result of decreasing T(p) with increas- 
ing temperature is a real effect. 

In the strong coupling phase with n/j/nj = 0.1, the 
trend of the spin relaxation is completely different. The 
spin relaxation T(p) is small and does not change much at 
low temperatures, then a peak structure appears just be- 
fore hitting the critical temperature. This non-monotonic 
trend can be explained by the energy shift of the Stoner 
continuum. In the strong coupling phase, the Stoner con- 
tinuum is separated from the gapless spin waves at low 
temperatures by a finite gap. The gap suppresses T(p) 
and makes it insensitive to temperature changes. The 
gap disappears at the temperature where the polariza- 
tion of the itinerant spins (a z ) is no longer fully polar- 
ized. Due to the intersection of the spin-wave dispersion 
and the Stoner continuum, enhanced spin relaxation is 
expected. An analysis of our numerical results reveal 
that the partial polarization of (a z ) and the peak in T(p) 
occur at the same temperature. 

One may notice the close relationship between the 
magnetization curve and the spin relaxation rate and 
their sensitive dependence on the density ratio nh/ni. 
This is reminiscent of the interesting differences between 
as-grown samples and annealed ones lj. The annealing 
process removes interstitial impurities and reduces the 
disorder strength. However, if the as-grown sample is al- 
ready in the diffusive regime, the annealed one can also 
be described by the virtual crystal approximation. We 
believe the major difference arises from the increase in 
carrier density and thus the density ratio n^/ni. The 
annealing process increases nh/ni from the strong cou- 
pling to the RKKY phase which changes the shape of the 
magnetization curve. Therefore, it would be exciting to 
study the spin dynamics for both as-grown and annealed 



samples and compare the experimental outcomes with 
our numerical predictions. Note that we did not include 
the realistic six-band Luttinger model here. In princi- 
ple, our approach can be generalized to include more 
bands by introducing more Green's functions but the self- 
consistent equations will be rather complicated, not fatal 
though. However, it is important to emphasize that the 
ratio rih/nj to enter the RKKY regime is expected to be 
larger, as compared to the estimate from the two-band 
model. 
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